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Abstract. 

The problem of Schrodinger propagation of a discontinuous wavefunction — 
diffraction in time — is studied under a new light. It is shown that the evolution map 
in phase space induces a set of affine transformations on discontinuous wavepackets, 
generating expansions similar to those of wavelet analysis. Such transformations are 
identified as the cause for the infinitesimal details in diffraction patterns. A simple 
case of an evolution map, such as SL{2) in a two-dimensional phase space, is shown 
to produce an infinite set of space-time trajectories of constant probability. The 
trajectories emerge from a breaking point of the initial wave. 
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1. Introduction 

What are the consequences of breaking a Schrodinger wave? In a seminal paper 
[1], Moshinsky studied the evolution problem of a quantum-mechanical wave blocked 
by a shutter and identified his results with Fresnel diflFraction. Since then, the so- 
called 'diffraction in time' has been widely studied [2j, and experimentally pursued [3j. 
Needless to say, the spatial side of the analogy (in the paraxial approximation [4J) has 
a long tradition of its own, mostly in the context of electromagnetic theory 

It is rarely recognized, however, that the discontinuities of the initial condition 
are the true origin of Schrodinger diflFraction as well as the peculiar intensity patterns 
accompanying the effect. The present paper delves into the subject by revealing its 
mathematical structure through the use of symmetries and self-similar relations obeyed 
by wave functions. It is shown that the afore mentioned patterns in space and time can 
be explained in a very general framework dealing with the evolution of discontinuous 
initial conditions. The role of symmetry in these quantum-mechanical problems will be 
crucial, as its infiuence in the evolution of the system shall give rise to non-classical 
trajectories of constant probability. 

Although diffr active effects are described by old and simple wave theories, the 
underlying explanation of the complexity in the emerging patterns is a subject of current 
discussion [7j and the resulting intricacies are the center of attention of other studies in 
connection with fractality [8j . It is worth to mention that the study of discontinuities in 
quantum-mechanical wave functions is also motivated by current physical applications 
ranging from molecular Talbot interferometers [9j to matter-wave optics based on cold 
atoms [TOlfn]. 

Structure of the paper: In the following section, we study the old problem of 
diffraction by edges, in order to show all the relevant features of patterns by means of 
elementary techniques, including self-similarity and caustics. In section 3 we present 
a generalization of the method and illustrate it with a few examples in section 4. 
These include the diffraction of a square packet in a parametric harmonic oscillator, 
the evolution of discontinuities under non-linear canonical evolution and their relation 
to wavelet expansions. As a final example, the evolution of a discontinuous packet 
under the Gross-Pitaevskii equation is shown numerically. In section 5 we give a brief 
conclusion. 

2. Diffraction by edges 

It is usually considered that, in the paraxial approximation, the problem of diffraction 
by edges in the x, z plane can be described by a plane wave in z > propagating parallel 
to the z axis, a discontinuous opaque screen placed at z = parallel to the X axis and the 
corresponding solution at z > 0, subject to the boundary conditions mentioned before. 
The solution at z > can be found by propagating the wave at z = 0, which is taken as 
the original plane wave but forced to vanish at the blocking screen i.e. a discontinuous 
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Figure 1. Probability density as a function of space and time. Two regions can be 
distiguished: 1) the near zone or short-time regime, where infinite oscihations and a 
main lobe can be identified, and 2) the far zone or long-time regime, where the wave 
function spreads 



initial condition. In the following we solve the corresponding Schrodinger propagation 
problem in natural units {h = 1 = m) using the time variable t instead of z. We have 



1, 



t) 



(1) 



where ipo has, in principle, an arbitrary number of discontinuities. See the single-slit 
example in figure [T] 

Let us consider a variant of the problem treated by Moshinsky [Ij, revisited by 
Nussenzveig [3j. Consider the intial condition 



^Jj{x,t = 0) = e{l/2- \x\) (2) 

corresponding to a square packet of unit width. The wave function at time t can be 
obtained via the free propagator fT2j , which is a Gaussian of the form 



The wavefunction reads 



^ix,t) 



-Erf 



1/2)/ V2t 



-Erf 



-i7r/4| 



x + l/2)/V2t 



(4) 



In this form we can obtain the intensity pattern in space and time, as shown in figure 
[TJ It is worth to mention that the properties of such a pattern can be inferred from 
the Cornu spiral, as it was pointed out long ago [Ij. However, the applicability of this 
geometrical picture is restricted to the appearence of Fresnel integrals in the solution. 
More general initial conditions demand the use of more powerful methods. 



Phase Space Evolution and Discontinuous Schrodinger Waves 



4 



2.1. Self- similarity in single- slit diffraction using special functions 

The origin of a self-similar pattern near the edges of the initial condition can be described 
alternatively by a replication formula , z.e. a relation describing the wave function in 
terms of itself. Using a trigonometric expansion of the initial condition 

1 2 1 

e{l/2 - \x\) = 2 + - E (7r(2n + l)(x + 1/2)) (5) 

and applying the Gaussian propagator to this expansion, we can prove in a very simple 
way the identity 



Erf 



6-^^^/2(1/2 - x)/^/t\ = J2 Cne'^-^''^'^ETi[r]n{x,t)] (6) 



n=0 



where we have used the definitions 



/c^ = 7r(2n + l), Cn = l/ikn (7) 
U^,t) = -klt/2 + kn{l/2 + x) (8) 
r)r,{x, t) = (1/2 -X- kj) /Vt. (9) 

Here we can infer the shape of the pattern near the edges by using a short time 
approximation of the functions on the r.h.s. of ([g]). We simply take the individual 
packets to be the initial conditions whose argument is displaced according to r]n{x^t)^ 
that is 



ij{x,t) =Y.Cn [e'^-^^''^e{l/2 - \vn{x,t)\) + e^«"(-^'*)^(l/2 - \Vn{-x,t)\) 

n=0 

+ ^^(1/2 -1x1) + 0(^2) (10) 

where 0{t'^) is a small correction for short times. It is interesting to note that in this 
expansion, the emergence of the pattern in space time can be seen as the superposition 
of many individual square packets moving along the trajectories 7y^(±x,t), with plane 
wave factors propagating along ^^(±x,t). The alternate signs indicate the contribution 
from each discontinuity or edge, while the index n parameterizes the velocity. This 
velocity also alternates sign according to the source of the rays, propagating to the left 
from the right edge and vice versa. The individual functions preserving their integrity 
on the r.h.s. of (10) can be regarded as building blocks, while the coefficients Cn weight 
the contribution of each block. For the resulting probability density plots, see figures 
[2| [3| [4] Furthermore, the building blocks can be replaced by any other localized shape 
such as that of a Gaussian or a tringular packet, and the resulting pattern generated 
by them will have the same features. See, in particular, figure [4j We distiguish the 
following properties: 
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• When the hmit of the sums is truncated (say n = rimax)^ the pattern has a finite 
number of osciUations near the discontinuities. 

• The second term Umax = 1 shows already a main peak: The original pattern to be 
replicated. Increasing n implies finer detail in the pattern near the edges. 

• Rays are described by (^^(±x,t) = const. The envelopes of such system of rays 
are caustics (Legendre transform) and can be obtained by eliminating n from the 
condition dn^ni^^t) = 0. It turns out that such envelopes are parabolas centered 
at the edges. See figure [5} 

• Classical methods from geometric optics produce a non-classical result. To this 
respect we should point out that the trajectories of (approximately) constant 
probability obtained from the caustics should not be interpreted as classical 
trajectories. In free evolution, the latter are simply rays, while the former are 
parabolas. 



The short time approximation used in (10) can always be improved. For example. 



in the case of initial space dependent phases, a velocity field given by the gradient 
of the phase produces additional motion on the initial packets or building blocks. 
Using the continuity equation for the initial probability density, one can infer the 
specific form of the motion and obtain the correct approximation. 
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Figure 2. Replication process using square packets as building blocks. The red circle 
indicates the pattern to be replicated. 
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Figure 3. Replication process using gaussian packets as building blocks. The red 
circle indicates the pattern to be replicated. 
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Figure 4. A comparison between the pattern obtained through the replication 
formula and the numerical calculation using Moshinsky functions. The expansions 
contain 30 terms. 



3. Generalizations 

Many of the results above can be produced in more general settings. Now that we have 
revealed the origin of the patterns, we may consider intial conditions with an arbitrary 
number of discontinuities, as they can be written again in terms of individual square 
pulses 



fill 
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Figure 5. The emergence of caustics as the culprits for the constant intensity hues in 
the pattern. As indicated in the text, the caustics can be computed exphcitly in the 
form of parabolas 

It is also possible to treat similar problems in many dimensions. Take, for example, the 
free evolution of a packet in three dimensions written as the product of pulses in each 
coordinate x^y^z. Even more general initial conditions with compact support can be 
considered, as long as we use the appropriate trigonometric expansions. 

In this section we discuss not only the generalizations of the method, but also the 
role of symmetry in the evolution of discontinuities. Most of our treatment rests on 
the fact that a wave function can be written in terms of itself and on the possibility of 
representing its evolution as the superposition of wave packets which move and leave 
their mark in space-time, but preserve their integrity in the process. 

3.1. Replication formula from canonical evolution 

A discontinuous wave can be written in terms of itself using step functions. We resort 
to a trivial but remarkable identity for an individual pulse: 



^o(^) = 0{l/2 - \x\)iJo{x). (12) 
The evolution can be computed straightforwardly 

ij{x, t) = UtM^) = Ute{i/2 - \x\)u.tUtMx) (13) 

= 9(1/2 -\x{t)mx,t) (14) 

where we have used the reversed Heisenberg picture to denote our evolved position 
operator x{t). As before, we use a trigonometric expansion of the step function in the 
form 

6(1/2 - \x(t)\) = t ^e±^'="(^«±V2) ^ y: (15) 
with a convenient abbreviation of the sum symbol given by 



/ oo 

n,± n=0 ± 



(16) 
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Finally, the wave can be written in terms of itself by noting that each exponential 
in the expansion is itself a unitary operator to be applied to the wave at the right. Such 
exponentials depend on the evolved position operator containing both x and p at t = 0. 
In this way, many coordinate transformations can be induced on ip: 

= i:^U^^M = 5:^e^«"(^'*)^(r)„(x,t),t) (17) 

There are many choices of x{t) for which induces coordinate transformations. As 
the most general way (in two-dimensional phase space) take x{t) as a function of x^p in 
the form 

x{t) = ft{x) + {gt{x),p}. (18) 

This operator induces a change of variables in x and an afhne transformation for 
each n in the trigonometric sums above. To fulfill the initial data problem, we set 
/o(x) — igQ{x) = X and go{x) = 1/2. Obviously, the function r]n{x^t) producing the 
change of variables is to be determined for our evolution of choice ( [Tsl ). Rotations, 
shearings and dilations can be obtained in this form. 

What we have achieved so far can be summarized as follows: 

• The elements Gt of an evolution group act on Hilbert space as Gf ; T~L — y T~L . 

• For ip G H with discontinuities in the form of pulses, one has ip = Qip. The 
operator is a projector. 

• In the reversed Heisenberg picture one has G^^^ = GtOG-t^ which admits an 
expansion in terms of unitary operators B^^^ = CnUn- 

• For each [4, a new transformation is produced: x^^^ = GtxG-t = ft{x) + {gt{x)^p}- 
The transformation is canonical, i.e.p^^^ = GtpG-t- 

• The overall result becomes ip{x^t) = ^nCn'^{Vn{x^t)^t). 

In the following we present some examples. 
4. Examples 

4^1. Diffraction in a parametric harmonic oscillator 

Take as elements of the evolution group, the linear symplectic matrices with blocks 
^4, C, D acting on phase space x,p In the simplest two-dimensional application, 
these blocks can be considered scalars with the property A{t)D{t) — B{t)C{t) = 1 and 
such that 



x{t) = A{t)x + B{t)p 
p{t) = C{t)x + D{t)p. 



(19) 
(20) 
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Figure 6. A comparison of diffraction patterns under a harmonic oscillator potential: 
Numerics, theory (with improved individual packets) and trajectories. 



After substitution and the use of Baker-Campbell-Hausdorff formula, the replication 
formula becomes 



t)= Y,-^ exp (^-iklA{t)B{t)/2 ± ikn{A{t)x ± 1/2 



X 



X ^{x±knB{t),t) (21) 

For short times, we have an expansion in terms of building blocks or a mother wave 
, given by the initial condition: 



—i / 

iIj{x, t)=Y.jr ± ikn{A{t)x ± 1/2 

X ^o{x±knB{t)) + o{t'^) (22) 

Furthermore, the short-time approximation can be improved by using a relation 
between harmonic and free evolution wavefunctions coming from the fact that both 
propagators are Gaussians. In the case of an initial condition given by the square 
packet, one obtains an interesting intensity pattern in agreement with the numerical 



evaluation. The trajectories are obtained from the phase factors in (22), parameterized 
by A{t)^B{t). For a time-independent harmonic oscillator, such functions are taken as 
A{t) = cos(cjt), S(t) = — 8m{(jjt) / (jj . See figure|6]for a comparison. 

4^2. The emergence of wavelets 

Another example of interest emerges from the application of squeezing operations on 
wave packets. This shall lead to expansions similar to those of wavelet analysis [14j. 
Consider the transformation 

x{t) = x+^-{x,p} (23) 
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p(t) = ^log(tp+l) (24) 

which is an unusual (but vahd) example of a canonical, non-linear transformation 
generated by a hamiltonian which can be written down explicitly: 

H = \{x,{l-e-'n/t'-p/t}. (25) 
The replication formula now reads 

(26) 

Which implies a true scaling of the wave packets by powers of the parameters. These are 
visible in the exponentials accompanying the argument x. Each power can be understood 
as a new scale factor 'zooming' into the pattern (this can be regarded as a microscope 
of increasing power). The function can be written for short times as 

^{x,t) = ^2 ^ exp (-Z7r(e^^-^ - l)x/t) e^^-^/^o (e^^^'/'x) 

(27) 

This is nothing else than a wavelet expansion, where the mother wavelet is given by 
the initial condition ipoi^)- 

4-3. Diffraction and the Gross- Pitaevskii equation 

Here, as a final example, we present the numerical solution to the evolution problem 
of a square packet governed by the Gross-Pit aevskii equation. 



-^V2 + ^|7/;(x,t)p 



V^(x,t) 



.dilj[x, t) 



'0(x,O) = '0o(x). 



(28) 



We have chosen distributions with slightly smooth edges, showing thus that the effect 
in question is robust. It is interesting to note that many of the features present in the 
linear Schrodinger equation, appear also in the non-linear case [15]. This occurs despite 
of the fact that in the numerical calculations, strong values of the coupling parameter 
g = 50, 100 have been used. See figure [7| 



A brief explanation of the effect in question can be given in terms of the Gross- 
Pitaevskii equation plus the previously obtained replication formula for free evolution. 
For 5^ = we have 



(29) 



n=0 
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Figure 7. Probability density for increasing values of ^ = 0, 50, 100. Courtesy of S. 
Arnold [15]. 

The non-linear potential can be incorporated in the evolution by means of the interaction 
picture. For short times and non-vanishing coupling, we have that the corresponding 
operator can be replaced by a simple phase factor in the propagation problem. This is 
indicated in the following expression. 

ijg{x, t) = r dx'K{x, x'] t, 0)e-^^'l^°(^)l Vo(^) + o{e) 

J —oo 

(30) 

which is in the form of Schrodinger propagation with an effective initial condition 
^-wA'^oix)? ijj^l^x) . The intitially space-dependent phases produce additional focusing or 
defocusing in the previously obtained patterns (linear case). In the case of positive 
coupling g we expect defocusing, preceded by a self-similar regime whose duration is 
shorter than in the free case. 

5. Concluding remarks 

In this contribution, devoted to the role of symmetries and self-similarity in diffraction 
patterns, we have shown that the Schrodinger evolution of discontinuities can be 
described by a wave function in terms of itself. This gave rise to a replicating pattern 
near the edges of an initial distribution. The notion of replication can be extended 
to other parabolic equations, including non-linear terms. A formal theory of the 
present treatment is to be developed: Envelope curves live in manifolds, rays in the 
exponents live in tangent spaces. It has been stressed that the symmetry groups 
producing the evolution in phase space induce transformations of wavepackets through 
the exponentiation of the evolved coordinate operators. 

It is desirable to find all these patterns in experimental setups, either using classical 
light, or the more ambitious manipulation of Bose-Einstein condensates. The current 
technology suggests [16j that the needed edges in the initial distributions can be 
produced under special configurations of the confining electromagnetic traps 
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